TDIR: Time-Delay Interferometric Ranging 
for Space-Borne Gravitational- Wave Detectors 



O 

o 

(N 
•*-> 

o 
O 
in 



> 

(N 



^t- : 
o ■ 

"o 
cr 

i ' 



X 



Massimo TintoQ Michele VallisneriQ and J. W. Armstrong 

Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109 

(Dated: February 7, 2008) 

Space-borne interferometric gravitational-wave detectors, sensitive in the low-frequency (mHz) 
band, will fly in the next decade. In these detectors, the spacecraft-to-spacecraft light-travel times 
will necessarily be unequal and time- varying, and (because of aberration) will have different values 
on up- and down-links. In such unequal-armlength interferometers, laser phase noise will be canceled 
by taking linear combinations of the laser- phase observables measured between pairs of spacecraft, 
appropriately time-shifted by the light propagation times along the corresponding arms. This pro- 
cedure, known as time-delay interferometry (TDI), requires an accurate knowledge of the light-time 
delays as functions of time. Here we propose a high-accuracy technique to estimate these time delays 
and study its use in the context of the Laser Interferometer Space Antenna (LISA) mission. We 
refer to this ranging technique, which relies on the TDI combinations themselves, as Time-Delay 
Interferometric Ranging (TDIR). For every TDI combination, we show that, by minimizing the rms 
power in that combination (averaged over integration times ~ 10 4 s) with respect to the time-delay 
parameters, we obtain estimates of the time delays accurate enough to cancel laser noise to a level 
well below the secondary noises. Thus TDIR allows the implementation of TDI without the use of 
dedicated inter-spacecraft ranging systems, with a potential simplification of the LISA design. In 
this paper we define the TDIR procedure formally, and we characterize its expected performance 
via simulations with the Synthetic LISA software package. 

PACS numbers: 04.80.Nn, 95.55.Ym, 07.60.Ly 



The Laser Interferometer Space Antenna (LISA) p| is 
a planned NASA-ESA mission to detect and study grav- 
itational waves (GWs) in the 10 _4 -1 Hz band, by ex- 
changing coherent laser beams between three widely sep- 
arated spacecraft. Each spacecraft contains two drag- free 
proof masses, which provide freely falling references for 
the GW-modulated relative spacecraft distances; these 
are measured by comparing the phases of the incoming 
and local lasers. It has been shown that the time series of 
the phase differences can be combined, with suitable time 
delays, to cancel the otherwise overwhelming laser phase 
noise, while preserving a GW response. This technique is 
known as time-delay interferometry (TDI; see 0, U and 
references therein). 

In the case of a stationary LISA spacecraft array, it was 
estimated Q that the time delays need to be known with 
an accuracy of about 100 ns, if the various TDI combina- 
tions are to work effectively, suppressing the residual laser 
phase fluctuations to a level below the secondary noises 
(such as the proof-mass and optical-path noises). For 
an array of spacecraft in relative motion along realistic 
Solar orbits, more complicated {second- generation) TDI 
combinations are needed; these require an even more ac- 
curate knowledge of the time delays [I| . The most direct 
implementation of TDI consists in triggering the phase 
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measurements at the correct delayed times (within the 
required accuracy), as suggested in Ref. Q. This ap- 
proach requires the real-time, onboard knowledge of the 
light-travel times between pairs of spacecraft, which de- 
termine the TDI time delays. Although the triggering 
approach is feasible in principle, it complicates the de- 
sign of the optical phasemeter system, and it requires 
an independent onboard ranging capability. Recently, it 
was pointed out [£j that the phase measurements at the 
specific times needed by the TDI algorithm can be com- 
puted in post-processing with the required accuracy, by 
the fractional-delay interpolation (FDI) 0,0 of regularly 
sampled data (with a sampling rate of 10 Hz for a GW 
measurement band extending to 1 Hz). 

In this communication, we show that FDI allows the 
implementation of a numerical variational procedure to 
determine the TDI time delays from the phase-difference 
measurements themselves, eliminating the need for an in- 
dependent onboard ranging capability. Since this varia- 
tional procedure relies on the TDI combinations, we refer 
to it as Time-Delay Interferometric Ranging (TDIR). 

Conventional spacecraft ranging is based on the mea- 
surement of either one-way or two-way delay times. In 
one-way ranging, two or more tones are coherently mod- 
ulated onto the transmitted carrier; the phases of these 
tones are measured at the receiver, differenced, and di- 
vided by the spanned bandwidth to yield the group delay 
and hence the time delay (up to an ambiguity of c divided 
by the spanned bandwidth of the ranging tones). In two- 
way ranging, a known ranging code is modulated on the 
transmitted carrier, which is transponded by a distant 
spacecraft back to the originator; the received signal is 
then cross-correlated with the ranging code to determine 
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the two-way time of flight. 

TDIR differs from these methods in that it uses the 
unmodulated laser noises in a three-element array, which 
are canceled in TDI combinations assembled with the 
correct inter-spacecraft light-travel times. This means 
that TDI can be used to estimate the light-travel times 



by minimizing the laser noise power in the TDI com- 
binations as a function of the postulated light-travel 
times: this process defines TDIR. As an example of how 
TDIR works, we shall here consider one of the second- 
generation TDI combinations, the unequal-armlength 
Michelson combination X\ Q, 



X, = 



(V31 + »7l3;2) + (V21 + ^12;3');2'2 + (V21 + ^12;3');33'2'2 + 0731 + ?7l 3 ;2);33'33'2'2 

(V21 + Vl2;3>) + (V31 + ?7i3;2);33' + fel + ?7i 3 ; 2 );2'233' + fal + »7l2;3' );2'22'233' 



(1) 



Here we use the notation of Ref. Q, where the j/y (for 
spacecraft indices i,j = 1,2,3, i ^ j) are linear com- 
binations of the inter-spacecraft phase measurements sy 
and of the inter-bench measurements ry made aboard 
the spacecraft, 



W21 = S21 — - [T32 — T\2 



V12 = S12 



Vl3 = S13 



V31 = s 31 + o [721 



[T32 — T\2\ 



[T21 ~ 731 



??32 = S32 



??23 = S 2 3 



731 



7"13 - T23J .T, 



7"13 — T '23j 



(2) 

and where indices prefixed by a semicolon delay observ- 
ables by the corresponding light-travel times L k , in the 
sequential order given (with the index k = 1 denoting 
the light travel time for the beam emitted from space- 
craft 2 and received at spacecraft 3, and likewise indices 
{2, 3, 1', 2', 3'} = {3 -> 1, 1 - 2, 3 -> 2, 1 -> 3, 2 -> 1}). 

The main contributions to the phase measurements sy 
and ry are given by 



s 31 — 013:2 ~ 031 + ^31 > 
S21 = 012:3' — 021 + S 2 1 
731 = 021 - 031 + A31 , 
721 = 031 - 021 + A 2 1 



(3) 



(and by cyclical permutations thereof), where the <^>y de- 
note the sum of laser phase fluctuations and of optical- 
bench motions (the former three to four orders of mag- 
nitude larger than the latter), and where the £y and 
the Ay denote the sum of all other fluctuations affecting 
the measurements, such as the secondary noise sources 
(proof mass and optical path) and GWs. 

The time-delay indices that appear in Eq. J3J represent 
the actual delays caused by the physical propagation of 
the laser signals across the LISA arms. By contrast, the 
hatted delays of Eq. need to be provided by the data 
analyst (or, in the triggering approach, by the onboard 
ranging subsystem) with the accuracy required for effec- 
tive laser noise cancellation. Thus, the Ai-based imple- 
mentation of TDIR works by minimizing the power in X\ 



with respect to the hatted delays L k . Since the TDI com- 
binations constructed with the actual delays cancel laser 
phase noise to a level 10 8 below the secondary noises 0, 
it follows that if we neglect all non-laser sources of phase 
noise (i.e., if we set £y = Ay = 0), the minimum of the 
power integral 



(4) 



will occur for L k — L k (with k — 1, 2, 3, 1', 2', 3'; here 
the superscript (°) denotes laser-noise-only quantities). 
The search for this minimum can be implemented in post 
processing, using FDI jy] to generate the needed sy and 
Ty samples at the delayed times corresponding to any 
choice of the L k . 

In reality, the presence of non-laser phase noises (pos- 
sibly including GWs) will displace the location of the 

minimum from L k uv:,i "" v ~~ v 



X[ n) (with 



X\' L ' obtained by setting all 
becomes 



Writing Xi = X? 

>ij = 0), the power integral 



I {n) (L k 



1 

T 



[Ai(L fc )] 2 ^, 



(5) 



or explicitly, 



(L /5 )=I<°>(L*) + ^y [X\ 



Mi 2 



dt 



T 



x[ n) xi \L k )dt. 



(6) 



(n) 

Here we have written the non-laser phase noise X{ as 
independent of the delays L k : this holds true for a search 
conducted sufficiently close to the true minimum, since 
the <pij are much larger than the £y and Ay, and so 
are their variations. The minimum of I^ n \L k ) can be 
displaced from L k — L k because the third term of Eq. (0 
[the cross-correlation integral of X± and x[°\L k )] can 
be negative and offset a concurrent increase in I^\L k ). 
The achievable time-delay accuracies will depend on the 
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level of the residual laser noise, the levels of the secondary 
noises in X± , and the integration time T. We expect the 
arm-length errors to be determined by the interplay of 
the first and third terms in Eq. ©. By equating the 
variance from the imperfect cancellation of the laser with 
the estimation-error variance of the cross-term in Eq. (fit, 
we can roughly estimate how well the time delays will 
be determined with TDIR: 8L k ~ (cr („) j 'a a/ p/T, 
where ov(«i and cr • ( i are the root-mean-squares of the 

secondary noises and of the time derivative of the laser 
noise in X\ , and p is the temporal width of the secondary- 
noise autocorrelation function. For nominal LISA noises 
and T ~ 10,000 s we thus expect 5L k of 30 ns or better 
to be achievable. 

An analogous technique was proposed by Giirsel and 
Tinto for the problem of determining the parame- 
ters of a GW burst observed in coincidence by a net- 
work of ground-based GW interferometers. In that case, 
a "phase-closure" condition was imposed on a family of 
linear combinations of the responses of three GW detec- 
tors. A GW burst would produce a zero response in the 
particular phase-closed combination corresponding to the 
source's position in the sky; the position could then be 
estimated by implementing a least-squares minimization 
procedure. In that case as well as for TDIR, the least- 
squares minimization procedure can be shown to be op- 
timal if the secondary noises have Gaussian distribution. 

We test TDIR for a realistic model of the LISA orbits 
and instruments by performing simulations with the Syn- 
thetic LISA software package [ljj . Because the present 
version of Synthetic LISA works in terms of frequencies 
rather than phases, we perform an analogue of the pro- 
cedure outlined above where all the phase variables are 
replaced by the corresponding fractional-frequency fluc- 
tuations. We generate a number of chunks of contiguous 
data for the Sij and measurements, sampled at in- 
tervals of 0.25 s, and containing pseudo-random laser, 
proof-mass, and optical-path noises at the nominal level 
set by the LISA pre-phase A specification [ij, [l(j. We 
consider chunk durations of 8,192, 16,384, and 32,768 s. 

The 18 noise processes (corresponding to the six lasers, 
proof masses, and optical paths) are assumed to be un- 
corrected, Gaussian, and stationary, with (respectively) 
white, /" 2 , and f 2 PSDs, band-limited at 1 Hz. The 
frequency-fluctuation measurements contain also the re- 
sponses due to GWs from two circular binaries with 
/gw — 1 and 3 mHz, located respectively at the vernal 
equinox and at ecliptic latitude 45° and longitude 120°. 
The strength of the two sources is adjusted to yield an 
optimal S/N of ~ 500 over a year (for X\), guaranteeing 
that there will be times of the year when each source will 
be clearly visible above the noise in an observation time 
- 10,000 s. 

We put the three LISA spacecraft on realistic trajec- 
tories, modeled as eccentric, inclined solar orbits with 
angular velocity ft — 27r/yr, average radius R/c ~ 499 s, 
and eccentricity e ~ 9.6 x 10 -3 The resulting time 
and direction dependence Q of the light travel times is 



then [HITlI 

Lk{t) =L+ i-(eL) sin(3fif - 3£o) 

- [^(ei) ± (QRL)] sin(fif - S k ) , 

where the plus (minus) refers to unprimed (primed) in- 
dices. In Eq. L/c ~ 16.68 s is the average light travel 
time, and 

4-7T 2tt 

(*i,<M3) = (£o,£o + -3-,£o + -3-), (8) 

with £o an arbitrary constant (set to in our simula- 
tions) giving the phase of the spacecraft motion around 
the guiding center of the LISA array. The starting times 
of the chunks are spread across a year to sample the time 
dependence of the L k and the directionality of the GW 
responses. 

Separately for each chunk, we minimize I^ n ' [L^it)] [Eq. 
starting from guesses for the Lk affected by errors > 
50 km/c, very much larger than typical accuracy of radio 
tracking from Earth [L|. The minimization is carried 
out using a Nelder-Mead simplex-based algorithm [l3j . 
The effective cancellation of laser noise with TDI requires 
modeling the time dependence of the travel times within 
the chunks. In our simulations, we use two such models: 

1. An orbital-dynamics model (ODM) given by Eq. 

(0, with eL, flRL, and £o taken as the indepen- 
dent search parameters with respect to which j' n J 
is minimized. We exclude L and f2 from the search 
because the dependence of the Lk(t) on such an ex- 
tended parameter set is degenerate on time-scales 
- 10,000 s. 

2. A linear model (LM) given by L k (t) =L Q k + L\(t- 
to) [with t set to the beginning of each chunk]. 
Because the expression for X\ does not contain the 
travel times L\ and Ly, our independent search 
parameters are the constants L\ and L\ for k = 
2, 2', 3, 3' (eight numbers altogether). 

Figures ^and [21 show the results of our simulations. The 
average travel-time errors AL displayed in Fig. are de- 
fined as AL = (AL 2 + AL 2 > + AL 3 + A£ 3 /)/4, with 




Because the noises have different realizations in each 
chunk and because the local behavior of the Lf-[t) [Eq. 
0] changes along the year, the average error AL of each 
chunk is a random variable. Its distribution is approxi- 
mated by the histograms of Fig. ^ which refer to popu- 
lations of respectively 512 (for T = 8,192 s), 256 (for T 
= 16,384 s), and 128 (for T = 32,768 s) chunks (hence 
the roughness of the curves) . 
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Distribution of light-time errors 







T = 8,192 S(LM) . 


0.20 


i ^ 


o\/n - 91 — i 

dvy — o.c. i 1 — ■ 






0.15 




T = 8,192 s(ODM) . 


0.10 




- 


0.05 
0.00 


^-i avg = 1 .82 


" — 






T= 16,384 s(ODM) . 


0.20 


^ avg = 1 .56 








0.15 




T = 32,768 s(ODM) . 


0.10 






0.05 
0.00 


avg = 1 .37 L_n_ 





1 2 3 4 5 6 

AL [m] 



FIG. 1: Distribution of errors AL [see Eq. JjJJ and the main 
text above it] in the determination of light travel times, using 
Xi-based TDIR with chunk durations of 8,192 s (for the LM 
and ODM models), and 16,384 and 32,768 s (for the ODM 
model only). As expected, the errors are lower for longer 
integration times T; for the LM model, the larger errors are 
due to the unmodeled curvature in the time dependence of 
the light-travel times. The distributions shown correspond to 
samples of 512, 256, and 128 chunks for T = 8,192, 16,384, 
and 32,768 s respectively, spread across a year. 



It turns out that the linear model is not quite sufficient 
to model the changes of the time-delays during the chunk 
lengths considered, since the minimum ALs [computed 
by least-squares fitting the parameters L° and Li to the 
L fc (t)] are in the range 0.25-2.60 m (for T = 8,192 s), 
1-10 m (for T = 16,384 s), and 4-40 m (for T = 32,768 
s). Thus, in Figs. ^ and [21 we show results only for the 
linear model with T = 8,192 s. [The minimization of 1^ 
over the LM parameters is delicate, because for X\ the 
laser- noise residuals turn out to depend strongly on AL2 , 
AL 3 / , and AL 2 < — AL 3 , but only weakly on AL 2 ' + AL3. 
In this case, the Nelder-Mead algorithm can be made to 
return accurate results by using the search parameters 

L®, L®,, L\, — L3, and L\, + L3, plus the corresponding 
L\ parameters.] 

Figure [21 shows the spectra of the residual laser noise 
[i.e., of x[ ^ at the minimum of I^(L),)], as compared 

with spectra of GWs and secondary noises [i.e., of X^]. 
The spectra are computed separately for each chunk us- 
ing triangle-windowed periodograms, and then averaged 
over the chunk populations. The two GW sources stand 
clearly above the secondary noises at 1 and 3 mHz. We 
see that the TDI cancellation of laser noise with TDIR- 
determined time-delays is essentially complete, with the 
residual laser noise several orders of magnitude below the 
secondary noises. We conclude that for T ~ 10,000 s, 
with the nominal LISA noises, and even in the presence 
of very strong GW signals, TDIR can easily reach the 
time-delay accuracy required for second-generation TDI. 



For frequencies below 10 mHz, the residual laser-noise 
power decays as / 6 , while the secondary noises decrease 
only as f 2 . We attribute the flattening near 0.1 mHz 
(which is insignificant with respect to the LISA perfor- 
mance) to a combination of leakage and aliasing in the 
numerical estimation of the spectra and of real effects due 
to the first non-constant terms in the travel time errors 
across the chunks. 

Finally, we estimate the power in the Fourier bins con- 
taining the simulated signals using two different time se- 
ries: in the first X\ was formed using perfectly known 
time delays, in the second using the TDIR-determincd 
time delays. Analyzing the 32,768-s chunks at the times 
along the simulated year where the signal amplitudes 
were maximum, we find that the signal powers in the 
two time series agree to the numerical precision of the 
calculation (about a part in 10 5 ). 

In summary, we propose a method that uses TDI 
and the intrinsic phase noise of the lasers in a three- 
element array to determine the inter-spacecraft light 
travel times. This method, Time-Delay Intcrfcrometric 
Ranging (TDIR), relies on the fact that TDI nulls all the 
laser noises when the time delays are chosen to match the 
travel times experienced by the laser beams as they prop- 
agate along the sides of the array. Simulations performed 
using the nominal LISA noises indicate that, for integra- 
tion times ~ 10,000 s, TDIR determines the time de- 
lays with accuracies sufficient to suppress the laser phase 
fluctuations to a level below the LISA secondary noises, 
while at the same time preserving GW signals. Our sim- 
ulations assume synchronized clocks aboard the space- 
craft, but we anticipate that TDIR may be extended to 
achieve synchronization, by minimizing noise power also 
with respect to clock parameters. 

TDIR has the potential of simplifying the LISA design, 
allowing the implementation of TDI without a separate 
inter-spacecraft ranging subsystem. At the very least, 
TDIR can supplement such a subsystem, allowing the 
synthesis of TDI combinations during ranging dropouts 
or glitches. TDIR may be applicable in other forthcoming 
space science missions that rely on spacecraft formation 
flying and on inter-spacecraft ranging measurements to 
achieve their science objectives. 

The TDIR technique presented in this paper was based 
on the TDI combination X\ . Since the accuracy in the es- 
timation of the time delays achievable by TDIR depends 
on the magnitudes of the secondary noises entering the 
specific TDI combination used, it is clear that it should 
be possible to optimize the effectiveness of TDIR over 
the space of TDI combinations Q. We will investigate 
this problem in a more extensive future article. 
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FIG. 2: Spectra of frequency laser noise (bottom curves) and of GW plus secondary noises (top curve) at the end of TDIR 
minimization using chunk durations of 8,192 s (for the LM and ODM models), and 16,384 and 32,768 s (for the ODM model 
only). We show averages of the spectra computed separately for each chunk using a triangle- windowed periodogram; the 
averages are taken over populations of of 512, 256, and 128 chunks for T — 8,192, 16,384, and 32,768 s respectively, spread 
across a year. In all cases, laser noise is suppressed to levels several orders of magnitude below the secondary noises: the cutout 
graph on the right shows that the typical laser-noise suppression factor with respect to secondary noise is ~ 5 x 10 3 for the 
worst case considered (8,192-s LM); it improves by a factor ~ 2 for 8,192-s ODM, and by factors of ~ v2 for each successive 
doubling of T. The GWs from the two circular binaries stand clearly above the noise at 1 and 3 mHz. 
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